Reveal non-Markovianity of open quantum systems via local operations 
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Non-Markovianity, as an important feature of general open quantum systems, is usually difficult 
to quantify with limited knowledge of how the plant that we are interested in interacts with its 
environment — the bath. It often happens that the reduced dynamics of the plant attached to a 
non-Markovian bath becomes indistinguishable from the one with a Markovian bath, if we left the 
entire system freely evolve. Here we show that non-Markovianity can be revealed via applying local 
unitary operations on the plant — they will influence the plant evolution at later times due to memory 
of the bath. This not only provides a new criterion for non-Markovianity, but also sheds light on 
protecting and recovering quantum coherence in non-Markovian systems, which will be useful for 
quantum- information processing. 



Introduction. — A systematic characterization of non- 
Markovianity in open quantum systems has recently be- 
come an active field of both theoretical and experimental 
study [IHH, which is mostly motivated by recent quests 
for robust quantum-information protocols. Understand- 
ing non-Markovianity can lead to novel strategies for pro- 
tecting the quantum coherent of the plant — the part of 
quantum system that we care about, e.g., the atomic 
spin. Especially when we have certain control over the 
plant, non-Markovianity of the bath, in principle, allows 
us to decouple the plant from the bath, which is known 
as dynamical decoupling ■ 

Although non-Markovian dynamics arises ubiquitously 
as long as the bath has a temporal response compara- 
ble to the dynamical time-scale of the plant, quantifying 
it systematically is not straightforward for general open 
quantum systems. In the literature, Breuer et al. pro- 
posed a measure based upon the evolution of the trace 
distance Tr|pi(i) — p2{i)\/2 between two different initial 
quantum states of the plant pi{Q) and p2(0) an in- 
crease of the trace distance gives a unequivocal signa- 
ture of non-Markovianity, as it indicates that information 
flows from the bath back to the plant. Alternatively, Ri- 
vas et al. proposed introducing an ancilla to entangle 
but not interact with the plant: an increase in the entan- 
glement between the plant and the ancilla during evo- 
lution signifies the existence of non-Markovian dynamics 
between the plant and the bath Q . These two important 
measures have been applied extensively in studying non- 
Markovian quantum systems, and has been compared 
theoretically [I^ and experimentally by Liu et al. 

We notice that the above-mentioned measures are fo- 
cused on the reduced dynamics of the plant which con- 
tains a limited amount of knowledge of the entire plant- 
bath system. It often happens that the law of time evo- 
lution of the plant becomes identical to that of a Marko- 
vian system, even though the bath still carries non-trivial 
information about the plant. Here, we propose a new cri- 
terion for non-Markovianity, assuming that we not only 
have access to the time evolution of the plant's density 
matrix, but can also carry out unitary operations on 
the plant (but not on the bath): the dynamics is non- 



Markovian if a local unitary operation on the plant at 
a given moment can influence the plant's law of evolu- 
tion at later times. We will use this criterion to reveal 
non-Markovianity in systems where the plant follows the 
time-local Markovian master equation: 

'pp{t) = -{l/h)[Hp{t),Pp{t)] + E^7^{t)^^Ppit) ^ (1) 

where 7i(i) > and Lindblad terms CiPp = 2AiPpAl — 
{AjAi, pp} with Ai being plant operators 11, 1^ . 

Reduced dynamics of the plant. — We first discuss some 
general features of the reduced dynamics of the plant, 
which help understand the new criterion. Suppose our 
plant-bath system evolves from t = Otot = T. We divide 
this process into N small segments with increment dt — 
T/N. The entire system undergoes a unitary evolution: 

PpbiT) = UN{dt) ■ ■ ■ Ui{dt)ppb{0)Ul{dt) ■ ■ ■ Ul,{dt) , (2) 

where ppb is the plant-bath density matrix and Ui{dt) is 
the unitary evolution operator of the total Hamiltonian 
at ti. The reduced dynamics for the plant is obtained by 
tracing over the bath at each step, as shown schematically 
in Fig.[l] and the density matrix of the plant evolves as: 



Pp{T) = ^N{dt) ■ ■ ■ <l,(di) • • • $i(dt)pp(0) , 



(3) 



where pp = TibiPpb] and the super-operator <I>,;(di) is a 
trace- preserving dynamical map at t ~ ti. 

In general, the dynamical map ^i{dt) relies on the his- 
tory of the plant-bath state Ppb{tk) {tk < ti)- For the 
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FIG. 1: (color online) A schematic showing how the reduced 
dynamics of the plant emerges from the full dynamics of the 
plant-bath system by tracing over the bath state at each step. 
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plant to be strictly Markovian, the bath's memory about 
the plant must not affect the plant's further evolution, so 
$i(dt) must depend only on the plant- bath state at ti. 
In the simplest case, <I>i(dt) is independent of time and 
the state — the dynamical map forms a semigroup with 
$(t + T) = $(t)$(T). The corresponding generator is the 
Lindblad super-operator C, namely $(t) = e^*, and the 
master equation for the plant is in the standard Lindblad 
form: Pp{t) = -{i/h)[Hp, pp{t)] + Cpp{t). 

In general, even though the dynamical map for 
a non-Markovian plant depends on the history of the 
plant-bath system before U , for most situations this time- 
nonlocal dependence can be consistently accounted for by 
a master equation for the plant that is local in time as 
shown by Chruscihski [l^ : 
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(4) 



Here the fact that Pp{t) is written in terms of Pp{t), 
i.e., the increment in pp only depends on pp locally, is 
merely a construction. The memory effect of the non- 
Markovian dynamics still persists in the sense that the 
form of dependence Ai actually varies from each initial 
plant-bath initial state ppb{0). It is in this context that 
the issue of "characterizing non-Markovianity" is raised, 
where one hopes the form of M provides a hint on the 
non-Markovianity of the system by which M is deduced. 
Unfortunately, as we shall show later, there exist situa- 
tions in which the form of Ai does not particularly differ 
from a Markovian master equation. 

Criterion for non-Markovianity. — A natural way to re- 
veal non-Markovianity is to explore the memory effect — 
the dependence of M. on /5pf,(0). This also has an oper- 
ational meaning: we should apply a time-local unitary 
operation on the plant-bath system and study whether 
and how it influences the plant evolution at later times. 
Let us restrict the operation to the plant (the bath is 
usually uncontrollable in actual systems). In terms of 
dynamical maps, our new criterion can be phrased as: 

The plant 's open quantum dynamics is non-Markovian 
if an instantaneous operation on the plant can lead to a 
change in its dynamical maps at later times. 

In other words, the plant dynamics is non-Markovian 
if Ppb{tk) Up(E)ib Ppb[tk)Ul ® h at tk < T, with Up a 
plant operator and 7^ identity for the bath, can lead to 

Pp{T) + $Ar(dO • • • 4,(d0 • • • $fc+l(d^)^7ppp(^fe)c>p^ (5) 

where $ are dynamical maps before applying Up ® 
This new criterion, to some extent, serves as an oper- 
ational definition for non-Markovianity. It explores the 
memory effect in the non-Markovian dynamics — in par- 
ticular, the dependence of the initial plant-bath state. 
As we will show in the following example, this criterion 
can indeed reveal non-Markovianity in systems of which 
the undisturbed evolution (before applying the unitary 
operation) is indistinguishable from the Markovian one. 




FIG. 2: (color online) A schematic showing the atom-cavity 
system. The cavity mode is coupled to the external continuum 
field (a Markovian bath) , and they together forms an effective 
non-Markovian bath for the atom. (Adapted from Ref. [l3|)- 



Example. — To illustrate this new criterion, here we 
consider the atom-cavity system as shown in Fig. [2] — a 
two-level atom coupled to a cavity mode which in turn 
couples to an external continuum — a quantum Wiener 
process that is equivalent to a zero-temperature Marko- 
vian bath |14l | . If we view the cavity mode and the exter- 
nal continuum together as the bath, the two-level atom — 
the plant — is effectively coupled to a damped cavity mode 
which is a non-Markovian dissipative bath, similar to the 
pseudo-mode model [isl . [lij . The corresponding Hamil- 
tonian for this system is given by [17| : 



ii =h{L0q/2)&z + ft Aa^a + hg{a^a) + a+a) 

+ ftV^[a6U0+a^&in(t)]. (6) 

Here is the Pauli matrix and a± = {ax^i(jy) j'i-] a and 
&in are the annihilation operators of the cavity mode and 
the input field of the external continuum with [a, a^] = I 
and [6in(i), b\^^{t')] = 5{t — t')] uiq is the atom transi- 
tion frequency and A is detune frequency of the cavity 
mode; g and 7 are the corresponding coupling constants. 
After tracing over the external continuum, the joint den- 
sity matrix of the atom and cavity satisfies the following 
Markovian master equation: 

p{t) ~ — i [{ujq/2)az + Aa^a + g{&-a) + a+a), p{t)\ 

+ -i[2ap{t)a} ~ {a)a, p{t)}]. (7) 

To further obtain the master equation for the atom by 
eliminating the cavity mode, we need to know initial state 
of the atom and the cavity mode. Under the usually- 
applied assumption, they initially are separable and the 
cavity mode is in the vacuum state: = IV'a) ® |0). 
As shown in Ref. [l7[ , the reduced density matrix of the 
atom pa{t) satisfies a time- local master equation: 

\ijJq 

Pa{t) - i \^^z +g3{/(t)}6-+(7_, Pa{t) 

+ .g5R{/(0}[2a_p„(t)a+-{a+a_, pa(t)}]. (8) 

Here the time-dependent function f(t) satisfies the fol- 
lowing Riccati equation: /(<) — i{yiq — A + ij)f{t) — 
9 P{t) = 9 with the initial condition /(O) = 0, and the 
solution is given by [l8l | : 
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(9) 
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+ z A)^ — ff^]^^^- In the tuned case 



with/3= [i(7- 
with Wg = A and a strong dissipation 7 > 2g, f{t) is real 
and positive, and we simply have: 



trace distance 



Pait) 



— f^z, Pa[t) 



gf{t)Cpa{t) (10) 



with Cpa{t) — 2(T_pa(t)(7+ — {(7+(T_, Pa{t)}. Such a mas- 
ter equation can also describe the case when the atom is 
directly coupled to the Markovian bath but with a time- 
dependent coupling rate, of which the Hamiltonian is: 

H = fiK/2)a, + h^Y{t)[a. bl{t) + a+ Si„(t)] (11) 

with 7'(t) = gf{t). Therefore, by looking at the unper- 
turbed evolution of the atom density matrix, we cannot 
tell whether the underlying dynamics is non-Markovian 
or not, even though the atom-cavity interaction is highly 
non-Markovian when 7 < Wg . 

If we perturb the atom by applying a unitary op- 
eration, its density matrix evolution will deviate from 
Eq. pUj) due to memory of the cavity mode. To see such 
a deviation, the most transparent way to look at the 
evolution of expectation values of the plant dynamical 
variables — (o'x), i^y), i^z), and to compare it with their 
Markovain evolution which is given by [from Eq. (|T0| ]: 



(12) 



Such a comparison is shown in Fig.|3] — the initial atom- 
cavity state is (|0) 4- \l))/\/2<» |0), Wg = A = .g = 1, 
7 = 2g, and a unitary operation on the atom: ct^ Ic 
is applied at t = I. As we can see, the Markovian and 
non-Markovian evolution are identical before the unitary 
operation, and deviate from each other after t — I. 

There are other figure of merits that can also quan- 
tify the deviation, e.g^ the trace distance between two 
different atom states |3| , and the entanglement strength 
between the atom and an introduced ancilla [3| . In the 
left panel of Fig. 21 we show the time evolution of the 
trace distance Tr\pi{t) — p2{t)\/2 with two initial states 
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FIG. 3: (color online) Plot showing the time evolution of 
{d'x.y,z}- The black curve shows the Markovian dynamics from 



Eq. (|12p : the red curve shows the non-Markovian dynamics 
derived from the total atom-cavity evolution given in Eq. ((Tjl . 
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FIG. 4: (color online) Plot showing the time evolution of 
the trace distance (left) and the concurrence (right). The red 
curve shows the non-Markovian evolution described by Eq. ^ 
and the black curve is the Markovian evolution from Eq. (|10|l . 



for the atom given by pi{0) = |1)(1| and /52(0) = 
In the right panel of Fig.|4l we show the time evolution 
of the concurrence — the ancilla entangled with the atom 
is another two-level system, and therefore the entangle- 
ment strength can be quantified by the concurrence [l9| . 
Initially, the atom and the ancilla are in the maximally- 
entangled state: ^[|0)|®|0) + |1)|® |1)]. Other specifica- 
tions are the same as those for producing Fig.|3l Both the 
trace distance and the entanglement strength increases 
after applying the local unitary operation on the atom. 
Such a revival of the quantum coherence clearly indicates 
non-Markovianity, just as expected for such a system. 

Dynamical recovering. — The increase of quantum co- 
herence just shown, when the plant is perturbed with 
local unitary operations, can be important for quantum- 
information processing. This allows us to recover infor- 
mation of the plant that is stored in the bath, which 
we can call "dynamical recovering" . More importantly, 
by combining local operations with the dynamical decou- 
pling protocols — applying a sequence of control pulses 0- 
ISl , we can characterize the plant-bath dynamics even for 
an unknown bath, which can help us find the optimal 
strategies for maintaining the quantum coherence. 

For illustration, we again use this atom-cavity system. 
In the left panel of Fig.O we show the time evolution of 
the trace distance for two different dynamical decoupling 
procedures: (i) the first one is a direct decoupling at t = 1 
by applying a sequence of az ® Ic within short intervals; 
(ii) the second one is a delayed decoupling (delayed by 
At = 0.1) after the local operation at t = 1. Their differ- 
ence in the trace distance tells us how much information 
of the atom is recovered from the cavity mode, and it 
varies depending on the delay time At. Interestingly, the 
maximal difference is achieved when the atom and the 
cavity mode becomes disentangled, as shown by the right 
panel of Fig. [5] — the concurrence for the atom-cavity en- 
tanglement vanishes when the trace distance achieves the 
local extremum poj . This is understandable — quantum 
information of the atom can on longer be recovered from 
the cavity mode once they are disentangled. 

For general plant-bath interactions, starting from any 
moment during the evolution, there should exist an op- 
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FIG. 5: (color online) The left panel shows the time evolution 
of the trace distance with a direct decoupling (black curve) 
and a delayed decoupling At = 0.1 (red curve) after applying 
(Tz 8 ic at f = 1; the right panel shows that the maximal 
recovering is achieved when the atom and the cavity mode 
are disentangled [red solid curve is the trace distance (same 
as the left panel) and the blue dashed curve is the concurrence 
for the atom-cavity entanglement]. 



timal sequence of local operations on the plant for maxi- 
mally recovering its information. Operationally, one can 
define a measure for non-Markovianity based on such dy- 
namical recovering. We introduce the set of plant state 
pairs with the initial trace distance D that is equal to 1: 
A = {{pi(0),/}2(0)}, Vi?{/3i(0),p2(0)} = 1}. Suppose at 
moment t, D\p\(t\pi{ty\ = a < 1, the measure is: 



Na ~ max 



1 — a 



(13) 



where l>c/,p(i) = TT^SU^ij) ® h ppb{t)U^{T) (g) h] is a 
sequence of unitary maps for the plant from t to t + t. 
Obviously, for Markovian systems, Na is equal to 0, and 
for non-Markovian systems in which the plant can recover 
its initial states via local operations, the measure is equal 
to 1. In general, Na ranges between and 1 depending 
on how strong the bath memory is, and also the moment 
t that we start to apply unitary operations. 

Discussion. — This criterion can be applied to study the 
recent experiment by Liu et al. In their setup, the po- 
larization degree of freedom of photons (the plant) cou- 
ples to the frequency degrees of freedom (the bath). Dur- 
ing the time evolution, these two degrees of freedom be- 
come entangled — the polarization undergoes decoherence 
when tracing over the frequency degrees of freedom. By 
changing the frequency profile of the bath, i.e., its quan- 
tum state, the time evolution of the trace distance for 
the polarization can either monotonically decrease or os- 
cillate. This was used to demonstrate switching between 
the Markovian and the non-Markovian regime. Interest- 
ingly, if a unitary operation is applied on the polarization 
by flipping it at a given moment, one will find that the 
quantum coherence of the polarization will revive at later 
times even in the so-called Markvoian regime, which in- 
dicates that the intrinsic dynamics is non-Markovian as 
the frequency degrees of freedom contains memory about 
the polarization. Therefore, this experiment can make a 
direct test of our new criterion. 



Conclusion. — We have presented a new criterion for 
non-Markovianity in general open quantum systems — the 
non-Markovianity manifests in terms of non-local change 
of the dynamical map after applying a local unitary op- 
eration on the plant. It allows us to tell whether a 
time-local positive map is only an artifact of a special 
initial quantum state of the plant-bath system or not. 
We have illustrated this criterion with the atom-cavity 
model. This work, on the one hand, helps clarify some 
subtleties of non-Markovianity in open quantum systems; 
on the other hand, it provides a route to probe the non- 
Markovian bath and to enhance quantum coherence in 
quantum-information processing. 
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